Anisotropic Domain Growth of ANNNI Model at Low Temperatures 
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We investigate the ordering kinetics for axial next nearest neighbor Ising (ANNNI) model in 
one and two dimensions by the multi-spin heat bath dynamical simulation. This dynamics enables 
us to overcome the pinning effect and to observe the dynamical scaling law for domain growth in 
the ANNNI model at zero temperature. The domain growth exponent is 1/2 isotropically both in 
the ferromagnetic and the dry- (commensurate) antiphase. In the wet- (commensurate) antiphase, 
however, it is approximately 1/3 in the modulated direction, whereas it remains 1/2 in the non- 
modulated direction. We suggest that these exponent values are dictated by 3 and 4 body diffusion- 
reaction processes of domain walls. 

PACS numbers: 05.70.Ln, 64.60. Cn, 02.70.Lq 



X 
S3 



The phase ordering kinetics of systems quenched from 
a high temperature disordered phase to a low tempera- 
ture ordered phase has been widely studied including the 
domain growth law and the dynamical scaling behavior of 
correlation functions It is generally accepted that 

the universal behavior of the domain growth depends on 
the presence of topological defects, the conservation law 
of the order parameter, and types of local energy bar- 
riers [^|-^. Virtually all the systems studied thus far 
including anisotropic models such as the axial next near- 
est neighbor Ising(ANNNI) model are believed to have 
self-similar domains and an isotropic domain growth law 
[0-|lO|. While previous studies of the ANNNI model did 
observe somewhat anisotropic dynamic structure factors 
at low temperature the conventional belief is that 

the exponents characterizing the growth are isotropic. 
The previous work used single spin update algorithms 
which cannot overcome the strong pinning effect at low 
temperature and leads to sluggish dynamical behavior. 
Also, in accord with the prior expectation of isotropic 
behavior, the total excess energy (in both direction) was 
used to find the growth exponent |^,^ . 

The principal theme of this paper is a careful numerical 
study of the dynamics of coarsening of the ANNNI model. 
Our results indicate isotropic growth at low temperature 
for a range of parameter values governing the competi- 
tion between the nearest and the next nearest neighbor 
exchange interaction in complete accord with previous 
expectations. However, we find in another part of pa- 
rameter space striking evidence for anisotropic domain 
growth of self-affine domains with the exponent values 
(1/2 and 1/3) along the different directions being con- 
trolled by diffusion-reaction processes of the domain walls 
(defects). These results are obtained using a new multi- 
spin heat bath dynamics(MHBD) algorithm which allows 
one to overcome the strong pinning effects even at zero 
temperature. Our work provides the first vivid demon- 
stration of anisotropic scaling behavior in the coarsening 
problem and establishes the power of the MHBD algo- 
rithm for effective equilibration of this rich system. 



The ANNNI model was first introduced to describe 
the equilibrium properties of spatially modulated struc- 
tures in magnetic and ferroelectric materials where the 
commensurate-incommensurate (C-INC) transition ex- 
ists [pd]-p^. The competing interactions in this model 
are ferromagnetic interactions ( Ji > 0) of nearest neigh- 
bor spins in all directions and antiferromagnetic inter- 
actions (J2 < 0) of next nearest neighbor spins in the 
modulation direction and give rise to a rich phase dia- 
gram: at low temperatures, depending on the competi- 
tion ratio K = —JijJx, the ferromagnetic phase exists for 
K < 0.5 and the (commensurate) modulated antiphase 
(• • • Ttiittii • • •) for K > 0.5. This antiphase consists 
of a wet (0.5 < k < 1) and a dry phase (k > 1), which 
are separated by a line of wetting transitions as the tem- 
perature increases. On introducing the effect of thermal 
entropy, the phase diagram shows a variety of structures 
such as an incommensurate phase and devil's staircases. 

Our MHBD algorithm considers a square block of 4 
spins in 2D or a string of 4 spins in ID, so that we update 
a spin-cluster conformation among the 16 possible states. 
Each state has its own probability in accordance with the 
Boltzmann weight in the heat bath algorithm. One of the 
16 states is selected randomly according to their probabil- 
ity of occurrence and consequently a new configuration 
of the 4 spins is obtained. We first applied MHBD to 
the simple Ising model on a square lattice of linear size 
L=:1000. Simulations were carried out at T = O.lJ/fcs 
up to 10^ Monte Carlo step (MCS), and averaged over 
50 initial configurations. The expected scaling collapse of 
the correlation function is obtained and is very consistent 
with that obtained with single spin update algorithms. 
Moreover, the analysis from the excess energy, the defect 
density, and the correlation function provides a domain 
growth exponent of 1/2 . Thus MHBD is confirmed to be 
an excellent approach for studying the ordering kinetics 
at low temperatures. 

We now consider a quenching of the linear ANNNI 
model, where the modulation of ( Ji > and J2 < ) 
interactions exists precisely along the axis , from high T 
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to T = because an ordered state exists only at T = 0. 
Right after a nucleation stage from the random initial 
states when k > 1, the system consists of !(• ■ • iti ' ' 

3(--aTTTi---), 4(---iTtTti---), 5(---iTTTTTi---), 

6(- • • iTTTTTTi ■ ■ O'^^^s of up or down spins in the 
sea of dimers. These 1, 3, 4, 5, 6-mers (called domain 
walls or defects) diffuse, annihilate and produce stable 
dimers by coUisons as time elapses. Because the den- 
sity of trimers would be higher than any other fc-mer 
except for the stable dimer, most of the excess energy 
of this nonequilibrium state comes from trimers. When 
three trimers merge together (• • • iTTTiiiTTTi ■ ■ ■) , 
a monomer is produced with 4 stable dimers around 
(• • • ITTiiTiiTTi • • •) and this monomer, when com- 
bined with a trimcr, will result in two dimers eventu- 
ally. This 3 body-collision and annihilation of trimers is 
the dominant decay processes of domain walls and pre- 
dicts the domain growth law L(t) ~ (i/logt)^/^, where 
L{t) is a characteristic length and the logarithmic cor- 
rection originates from the fact that D=l is the criti- 
cal dimension for the 3 body diffusion-reaction process 
1^,^. However, for 0.5 < k < 1, the domain walls 
are strongly pinned if their motion is controlled by a sin- 
gle spin update and thus such 3 body annihilation or 
other decay processes of domain walls can never occur. 
But, applying MHBD with a string of 2, 3 or 4 spins, 
4 body annihilation process of trimers can occur which 
then becomes the dominant decay process in this regime: 
(3333) becomes (2433), (2442), or (2622) with the same 
energy, and (2622) becomes 6 dimers with an associated 
energy decrease AE = 4Ji + 8J2 < for 0.5 < k < 1. 
So, the rate equation for the density A{t) of trimers, 
^^^^ = —4:A{t)'^, is obtained which immediately predicts 
the growth law L{t) - A{t)-^ - t^/^ 

When 3 or 4 body collision of trimers occur, 4-mer, 5- 
mer, and 6-mer can be considered as the resonance of two 
trimers, three trimers, and four trimers. When 3 body 
coUison of trimers is the dominant decay process among 
that of fc-mer domain walls, the density of fc-mer domain 
walls will have a descending order by 3, 1, 4, 5, 6-mers, 
whereas when 4 body coUison of trimers is the dominant 
decay process, it will have a descending order by 3, 4, 5, 
6, 1-mers. Therefore, it is advantageous and important 
to look at the most and the next most dominant decay 
process of domain walls in the system, which predicts the 
domain growth law of either i^/^ or t^^^ directly. 

We simulated a phase ordering of ID-ANNNI model 
employing MHBD with blocks of either 2 or 4 spins for 
a system size 4000 up to lO'^ MCS at k=0.2, 0.4, 0.49, 
0.51, 0.6, 0.9, 1.0, 1.01, 1.1, 1.5, 2.0, 3.0, 5.0, and 10.0. 
The densities of fc-mer domain walls, the excess energy 
as well as the correlation function are calculated at each 
time and all quantities are averaged over an ensemble of 
2000 runs. Figure 1 shows the densities of fc-mer domain 
walls at K — 0.9 and k — 1.1 in the (commensurate) 
antiphase. In the wet phase {k — 0.9), the descending 
order in their densities are 3, 4, 5, 6, 1-mers and in the 



dry phase {k = 1.1), they are 3, 1, 4, 5, 6-mers. There- 
fore, the former and the latter ought to have the domain 
growth law of t^/"^ and t^/^, respectively. These results 
are self-consistently supported by our analysis of the two- 
point correlation function. 
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FIG. 1. The densities Ai^(t} of fc-mer domain walls for the 

ID-ANNNI model (a) in the wet phase (k = 0.9) ; 3,4,5,6,1-mer 

from top and (b) in the dry phase (k = 1.1) : 3, 1,4, 5, 6-mer from 

top. 
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FIG. 2. The scaled correlation functions of ID-ANNNI model in 
the ferromagnetic(ft = 0.4), the dry(K = 1.1), and the wet (k = 0.6) 
phase from top, which are collapsed for different times t = 80, 160, 
320, 640, 1280, 2560, 5120, 10240 MCS at each k values. Inset 
shows the growth exponents in each phases. 

The two point correlation function C (?', t) is defined as 
C(r, t) =< Si ■ Si+r ■ (—1)''/^ > suitable for the antiphase 
, where r is even number, and C(r, t) =< St ■ Si+r > for 
the ferromagnetic phase. The growth exponent n is eval- 
uated by the decay of the excess energy or the defect den- 
sity and by the scaling collapse of the correlation function 
C{r/L{t)), where L{t) is obtained from C{L{t),t) = 1/2 
(see Fig. 2). For the ferromagnetic phase (k < 0.5), 
we obtain n = 0.5 by fitting to a power law L{t) ^ i". 
For the dry-antiphase (k > 1), n ~ 0.46 on fitting to 
L{t) ~ while n ~ 0.52 on fitting to L{t) ~ {t/\ogt)"-. 
In order to detect the logarithmic correction, more exten- 
sive simulations with system size 2 x 10* were performed 
up to 6.5 X lO'^MCS at K — 0.9 and 1.1 averaged over 
300 samples. The log-log plot of L{t) versus t at k = 1.1 
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is not a straight line at longer times but has an upward 
curvature. On the other hand, the log-log plot of L{t) 
versus t/\og{t) shows a nice straight line. This is fully 
consistent with the fact that D=l is the critical dimen- 
sion of 3 body diffusion-reaction processes For 
the wet-antiphase (0.5 < k < 1), the growth exponent 
n ~ 0.35 with a power law fitting, which is close to 1/3 
and consistent with the theoretical prediction . This 
is a rare realization of 4 body diffusion-reaction process 
whose dynamic scaling is first observed in our numerical 
simulations. 

For a 2D-ANNNI model, the ground state in the an- 
tiphase has 4 degenerate states (A,B,C,D). If the ground 
state at two boundaries along the modulated (y) direc- 
tion of the finite system is fixed by A-state at one end 
and D-state at the other end, B and C-states can exist in 
between at equilibrium in the wet phase (0.5 < k < 1). 
In the dry phase (k > 1), the A-state meets the D-state 
directly without the appearance of B and C-states. The 
domain patterns in the dry and the wet phases are quite 
different (see Fig. 3): For a dry phase, it looks like a typ- 
ical domain growth pattern having 4 degenerate ground 
states. But, for a wet phase, it is elongated (growing 
faster) along the non-modulated {x) direction. At T = 
the wetting transition, similar to that of liquid between 
air and substrate, occurs at k — 1.0, and there exists 
a line of wetting transition as the temperature increases 
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3. The domain patterns of 2D-ANNNI model after 



0.6 (the wet phase) (b) at /t = 2.0 (the dry 



We performed MHBD simulations for ordering kinet- 
ics of the 2D-ANNNI model at T = for several n 
values (k=0.2, 0.4, 0.6, 0.9, 1.0, 1.1, and 2.0) taking 4 
spins as a square block. The system size is 2048 x 2048 
and 30 samples are accumulated up to lO^MCS. MHBD 
could remove the pinning effect for k < 1 so that do- 
mains grow well in both the ferromagnetic and the an- 
tiphase. The two point correlation functions for the an- 
tiphase were defined separately by Cj;(x, i) = < Si-Si+x > 
in the non-modulated {x) direction and Cy{y,t) —< 
Si ■ Si+y ■ (—1)^/^ > in the modulated (y) direction. The 
excess energy and the length of domain walls (defects) 



are calculated separately for each direction. The densi- 
ties of /c-mers along the y-direction are also calculated as 
in the ID case. The growth exponents are evaluated by 
Lx{t) ^ t""" and Ly{t) ^ obtained from the scaling 
collapse of correlation functions and by the domain wall 
densities Px{t) ^ t~^'"' and Py{t) ^ 

The growth exponents Ux and Uy are 1/2 at k = 0.2 
and 0.4 (ferromagnetic phase) and also at k = 1.1 and 
2.0 (dry-antiphase). But, at k = 0.6,0.9,1.0 (wet- 
antiphase), Ux and Uy are anisotropic: Ux ~ 0.55 and 
Uy ~ 0.39. One may wonder whether these values are 
transient and may both approach 1/2 (or 1/3 and 1/2 re- 
spectively) in the very long time limit. More simulations 
with the larger system size 4096 x 4096 were performed up 
to 2 X lO'^MCS at K = 0.8 and 1.0. Further more, MHBD 
with 3x3 and 4x4 block cells were used under the same 
conditions. However, the exponents Hx — 0.55 and Uy ~ 
0.39 do not change significantly ||l^. Hence one may also 
expect that the correction term must be included in Lx (t) 
and Ly(t). We considered two kinds of corrections: The 



first is p{t) 



L{t) ~ t"(l + t-^) and the second is 



~ t"(logi)". The non-linear fitting of the former does 
not give consistent results for different time intervals and 
K values. However, the same fitting of the latter gives 
Ux — 1/2, a ~ 1/2 and Uy ~ 1/3, a ~ 1/3 for many dif- 
ferent time intervals and k values. So, the anisotropic 
domain growth laws in the wet-antiphase are consistent 
with Lxit) ~ (i log t) and Ly{t) ~ (tlogt)!/^. Fig- 
ure 4 shows the nice scaling collapse of correlation func- 
tions both for the non-modulated (a;) direction and for 
the modulated (y) direction in the wet and dry phase. A 
similar collapse is observed for a plot of the domain wall 
density (not shown). 



o 

c 

CO 




O 

"O 

c 
ns 




x/L,(t)ory/L (t) x/L,(t) or y/L (t) 

(a) (b) 
FIG. 4. The scaled correlation function of the 2D-ANNNI 
model for non-modulated (x) direction ( top curve : C{x/ Lx{t)) 
vs x/Lx{t) ) and for modulated (y) direction ( bottom curve : 
C{y / Ly{t)) vs y/Ly(t) ) (a) in the wet phase (k=1.0) and (b) in 
the dry phase (k=2.0) collapsed for different times t = 40, 80, 160, 
320, 640, 1280, 2560, 5120, 10240 MCS at each re values. The 
similar collapse of correlation functions in the ferromagnetic phase 
was also achieved. Inset shows the growth exponents in different 
phases. 



The result from the densities of fc-mers along the y- 
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direction further supports and is consistent with the 
above remarkable result because their descending se- 
quence is 3,4,5,6,1 at K = 0.9 and 3,1,4,5,6 at k = 1.1 
in agreement with the ID case (see Fig. 5). So, one 
may identify that the 4 (3) body diffusion-reaction pro- 
cess of domain walls is important, as in ID, as the domi- 
nant decay process in the wet (dry) phase of 2D-ANNNI 
model, which results in the dominant domain growth law 

of (^1/2). 



phenomena in the (3 x 1) phase. It will be worthwhile to 
probe a temporal evolution of structual factor for these 
systems experimentally, and compare with our results. 
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FIG. 5. The densities of fc-mer for the 2D-ANNNI modeL (a) 
K = 0.9 (the wet phase) : 3,4,5,6,1-mer from top. (b) ft = 1.1 (the 
dry phase) : 3,1,4,5,6-mer from top. 

One may actually consider that the domain growth at 
T = is possible since MHBD is adopted. However, 
the simulation by a single spin update, Glauber dynam- 
ics, at T = OAJi/ks in the wet phase gives the similar 
anisotropic growth exponents Ux — 0.56 and Uy ~ 0.40. 
The two point correlation function along the modulated 
(y) direction has also the same dip as that from using 
MHBD. Therefore, a single spin update dynamics can 
give the similar but less precise result for T > 0, since it 
severly suffers from the strong pinning effect at the low 
temperatures. 

We have presented a summary of our results of de- 
tailed investigation of the ordering kinetics of the ID 
and the 2D-ANNNI model at T = using the multi- 
spin heat bath dynamical simulation to overcome the 
pinning effect. Our simulations show that the dominant 
domain growth exponent becomes 1/2 isotropically both 
in the ferromagnetic and dry- (commensurate) antiphase. 
In the wet- (commensurate) antiphase, however, it is ap- 
proximately 1/3 for the modulated direction, whereas it 
remains 1/2 for the non-modulated direction. The ex- 
ponent values are explained by 3 and 4 body diffusion- 
reaction processes of domain walls (defects) . Contrary to 
conventional belief regarding universality, our data sug- 
gest an extremely unusual situation in which the very na- 
ture of the domain growth law critically depends on the 
ratio between the strength of competing interactions. 

In the context of domain growth dynamics of exper- 
imentally relevant systems, the chemisorbed adsorbate 
systems such as O/Pd(110) and H/Fe(110) are of partic- 
ular interest |l6| , [l9| , p0| . Both systems exhibit a wetting 
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